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Abstract 



Bound states, such as qq and qq, may exist in the Quark Gluon Plasma. As the system is at high 
density, the volume of the bound states may evoke a reduction to the phase space. We introduce an 
extended bag model to investigate qualitatively the volume effect on the properties of the system. 
We find a limit temperature where the bound states start to be completely melted. 

PACS number(s): 12.38.Mh 25.75.-q 
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1. Introduction 

It is discussed recently that the Quark Gluon Plasma (QGP) may contain bound states, especially 
binary bound states such as qq, qq, gq and gg at temperatures above the critical value of deconfinement 
phase transition Tc CO EI El ■ Inserting these new degrees of freedom into the QGP will change its properties 
such as pressure and energy density. From lattice calculations jj, the masses of quarks and gluons are 
very heavy, such as niq^g/T '^SatTftil.S — 2Tc. If the bound states with finite volumes are introduced 
to QGP as well as the point like quarks and gluons with current masses, the high pressure at high T will 
be easy to be reached. 

The influence of particle volume is widely discussed in the Van der Waals Excluded Volume Model 
(VdW) to describe hadronic gas at high density. In QGP, the number density is very high so that 

the ball-like bound states would be quite close to each other as discussed by J. W. Clark, J. Cleymans 
and J. Rafelski [?]■ The volume effect of bound states may change the properties of such strongly coupled 
system remarkably. 

In this paper, detailed calculations will be presented to describe the influences of binary bound states 
in QGP. This may be helpful for us to understand the Equation of State (EOS) of QGP, which is quite 
important for the study of relativistic heavy ion collisions^. Many-body bound states such as baryons 
may also exist P, but they are neglected here for their large masses. 

2. Influence of the Bound State Volume 

Generally, statistical values such as e and P could be described as corresponding integrals, 



where fi{p) is the momentum distribution of number density for i particles. If all the components are 
point particles, the integral of coordinates equals to the system volume. When particles with excluded 
volume exist in the system, some of the coordinate space will be occupied. 
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Regarding those particles as rigid balls, the integral can described as, 

rf'r = (1 - yop)V, (2) 



where Tjop is the ratio of the volume occupation of the components. Then, for the bound states with 
volumes, one has 



Fi = (1 - Vop)FiQ 

- i^-yop)-^ lYMp)Mp)d'p, (3) 

This situation is very similar to the equation of Van der Waals in terms of the influence of particle 
volumes, {V — Nb), in 



{V - Nb) = NkT. (4) 



In relativistic conditions, j/op is defined as, 

Vop = (1 - Vop)yopO- (5) 

Then, 

Vop = T^^, (6) 
1 + 2/opO 

where [TO]. 

Here, Vi is the rest volume of particle i depending on temperature or energy density. ~ 7j^^ is the 
factor of Lorenz-contract. 

For point particles such as quarks and gluons, the "rigid balls" are transparent in fact. That means 
quarks and gluons can pass through the balls, but not completely freely as in vacuum. This situation is 
equivalent to a reduction of phase space, especially in low momentum regions. 

As an assumption, a parameter Cq. < 1 is used to describe the transparency for all the balls. Cq. is 
supposed to be a function of T or other variables, but the details are not clear. Approximately it would 
remain constant at given temperatures and is regarded as an adjustable parameter in our calculations to 
provide a preliminary demonstration. As a comment, the scattering of quarks and gluons with the bound 
state, which would destroy some of the bound states, will occur in the same time. This effect is part of 
the chemical equilibrium in QGP. 

Then, for point particles, 

F^^[{l-Vop) + Co,Vop]F,o■ (8) 

Specially, P20 a-nd £20 are used to describe the pressure and energy density of binary bound states 
and Pio, £10 are used to describe those of quarks and gluons. Then for the whole system, 

f = (1 - yop)P20 + [ (1 - Vop) + CaVoplPlO, 

£ = - y op)£2o + [ -y op) + Cayop]£ 10- (9) 

In order to study the influence on Quark Gluon Plasma when bound states with volumes are con- 
sidered, the expressions of masses and volumes of those binary bound states and their dependence on 
T or £ are needed. Generally, bound state masses would increase with temperature or energy density. 
This property makes them fewer and fewer, less stable and less important in QGP at high temperatures 
[3]. Volumes are also supposed to increase at the same time. A phenomenological model is required to 
construct masses and volumes in such styles. One of the suitable models to build volume is the Bag 
Model^IJIZj. It will be extended approximately in the condition of finite temperature. 
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3. Extended Bag Model 

MIT Bag Model[TT] describes the quark confinements in hadrons by additional boundary conditions. 
The mass of a hadron is expressed as ^21 > 

M{R)^Ev + Eo + Eq + Em+Ee, (10) 

where Ey = BVb is the volume energy, Eq = —Zq/R is the zero-point energy with constant Zq and Eq 
is the contribution of the confined quarks in the spherical region. Eq will decrease to the sum of quark 
masses when the radius of the bag expands. Em and Ee are the color magnetic exchange energy and 
the color electric interaction energy. By the quadratic boundary condition of pressure equilibrium, which 
is equivalence to minimizing M{R) with respect to R. Then, the bag radius i?o can be determined. Bag 
parameter was fit to be B^^^ — 14:5MeV12, . 

Bag model is originally used to describe the free hadrons in vacuum. When dense medium at high 
temperature is considered, the model need to be extended. Some discussions based on GCM or NJL 
model in conditions of finite density |15| noted that B would decrease with temperature or medium 

density. 

On the other hand, when discussing the phase transition between hadron matter and QGP, B is also 
introduced to balance the pressure difference between different phases. The value is fit to be i?^/** = 
230MeF2ni at Tc = WAMeV, and there is an approximate relation of i? oc Tc"*. 

These two directions can hardly meet together. It is also noted that B is originally defined to be both 
energy density difference and pressure difference between the opposite sides of the bag boundary. When 
medium exists, they will not be the same any more. The linear boundary conditions and the quadratic 
boundary conditions|ll) will also change. It means the model requires a small extension. 

In the extension, B = Bq is only regarded as the contributions of the sea quarks and gluons in basic 
states, which on whole contain zero quantum numbers. Bq is same to that in the original MIT Bag Model 
in vacuum. Keeping the pressure equilibrium or distinguishing different vacuum types between both sides 
of bag boundary by B of the same value may not be needed, as the color confinement may not exist in 
QGP. After all, this extension is just used to carry out a mechanism to discuss volumes and masses. The 
minimizing condition now is used on the whole energy of the system including medium instead, 

E{R) = Esag + Emed 

= BqVb + Ef) + Eq + Em + Ee + emedout{Vo — Vb) + SmedinVB, 
— [^0 — {^medout ~ £medin)]VB + Eq + Eq + Em + Ee + SmedoutVo, 

where Vq is the system volume and (emedotitVo) is constant to R. At the same time, the effective bag 
parameter B^ff — Bq — {Smedout — Smedm) decreases, which is consistent with the opinions of 13 531 CHI- 
While, the bag mass is still defined as Eq IIUII with Ey = BqVb- Direct dependence on chemical potentials 
is neglected. 

However, it is noted that Smedout and Smedin are not the summation of the energy per unit volume 
for the particles located inside or outside the bag boundary. That is because the volume effect can only 
restrict the coordinate space of the particle barycenters, instead of the particle wave functions or mass 
distributions. So, the particles located in one side of the boundary could contribute the energy density 
appropriately to the other side. Qualitatively, 

barycenter barycenter 

in out 
barycenter barycenter 

in out 

where, and e^^ are the energy per unit volume for the particles located inside and outside the bag 
boundary, is the effective wave function. Considering that the energy density difference between 



1 

^medin — 

^in 



1 

^medout — TT 

^ nil. 
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boundaries would be much smaller than expected by the superimposition of particle wave functions, a 
rest factor Cij could be introduced as a kind of average for a specific bag j, 



Summation to i and neglecting the difference of j, one has approximately, 

{^medout ^raedin) j ~ ^^med^j ~ ^^^med ~ /(^^med)i (-^-^) 

where, 

A£,„ed = e^^ - = (1 - yop)[e2a + (1 - Ca)£io], (12) 

f{Aemed) should increase slowly with Asmed or T, but the detailed expression may be less important 
for the final qualitative results like the Ising model. To avoid the divergence of numerical calculations, 

cAe^e. - Bo (^^1^) « Boil - e'T) = Bo{l - e-<^=^-"-). 

After defining, 

fie) ^ Boil (13) 
the bag masses and volumes could be determined by minimizing 

EiR) = (Bo - fiAe^ed))VB + Eo + Eg + Em + Ee- (14) 



- S Diquark 
-A Diquark 




Figure 1: Diquark masses with two light 
quarks via /(Aemed), Bo in unit. The up- 
per curve is the mass of axial-vector diquark. 
The lower curve is the mass of scalar diquark. 
Where, B^^'^ = U5MeV, Zo = 1.84 and = 
0.55 as unadjustable constants from |12| . 



Figure 2: Masses of color singlets ("mesons") 
with a light quark and a light anti-quark via 
/(Acmed), Bo in unit. The upper curve is the 
mass of vector bound state. The lower curve is 
the mass of pseudo-scalar bound state. Where, 
Bo^'^ = UbMeV, Zo = 1.84 and = 0.55. 



Masses of some diquarks (gq)|E| and qq singlets are shown in Fig |^ and They increase slowly 
with fiAsmed), until it is close to Bo, which makes the bag masses expand to infinite. The bag volumes 
also increase as the same styles. It is noted that in QGP with bound states, where the color and chiral 
symmetry is completely or partly restored, the pseudo-scalar qq singlets may not be Goldstone Bosons. 
They are considered much heavier than the pions in vacuum. 
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-In( 1 -f(Ae^ J/B, ) Ae'(e,„(1 -Ca)(1 +300/T)) 



Figure 3: Masses of pseudo-scalar singlets 
("mesons") and scalar diquarks with different 
components via f{Asmed)- When /(Ae^erf) is 
high enough, bound states with heavier quarks 
become smaller than those with light quarks. 



Figure 4: Demonstration of the solutions of 
the self-consistent equation where Ca = 

0.8 and Ce = 5 X IQ-^^MeV-^. Factor eio(l - 
C'a)(l + 300/r) is just used for the numerical 
calculations, where eio is the original energy 
density of quarks and gluons without bound 
states, see in Eq 



Although masses of bound states with light quarks are smaller than those with heavier ones in vacuum, 
they will be larger when f{A£med) grows high enough close to Bq. Demonstrations oi S ~ bound states 
are shown in Fig ||2Jl, 5 = 1 bound states act as the same styles too. This implies that ss and ss will 
be more important and more stable in conditions of high temperature in QGP relative to those bound 
states with light quarks. It is also consistent with the opinions that heavy quark bound states, such as 
J/vE", will contain larger effective potentials and gain more stabilities in QGP, although charm quark can 
not be used in this model so far. 

4. Statistical Properties of QGP with Bound States 

There are 394 binary bound states all together for Nf — 2 system and 588 for Nf = 3 system. In 
our calculations, it is supposed there is no color SU(3) symmetry breaking in Quark Gluon Plasma. 
That means some antisymmetric states with two identical particles such as J — 1 gg states and some of 
pseudo-scalar diquarks do not take into account. 

The system is regarded as an "infinitive" system approximately, with low chemical potentials. But 
leptons, which will escape from the finite region created by the relativistic heavy ion collisions, are 
neglected. 

Two light quarks and a heavy strange quark is used in the discussions, 

rUg = OMeV, mu,d = OMeV, = llOMeV] 

The current masses are supposed to remain still in the calculations when the temperature changes. 
Respectively, the components are considered as ideal Fermi or Boson gases, 

hip) = E^Zl ' 

e T ±1 

Combining Equations © © Q) © (^Ul) OS OS and strange quantum number conservation, one gets 

rrii — rn(A£), 

Vr = F(A£), 
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Hop = yop{mi,Vi,T,fiB,f^s), (15) 

Us = iis{mi,yop,T,HB,tJ^s), 
where Ae = lS.emed a-nd /is is the strange quark chemical potential, 

Eqs p5|l is a set of self-consistent equations equivalent to the forms with /i^ and Ae below, 

/is = <^i(T, Ae, /iB,Ms): 
Ae = </32(T, Ae, /iB,/is), 

Setting /is = 0, one has, 

Ae = ^(T,Ae), (16) 

In the calculations, gg states are neglected as their heavy masses as " Glueballs" |18|. Contributions of qg 
states are deduced by treating them as specific diquarks. For simplicity, Em and Ee are neglected from 
Eq (|14|l . when particle radii are larger than 16 GeT^~^. 

It is noted that there are only two adjustable parameters, Cq and Ce, in our calculations, which makes 
the results more comprehensive. A demonstration of the solutions of Eq (|16|l is shown in Fig The 
ratio of occupations and energy density from the solutions under different parameters are presented in 
Fig (0)-®- As shown in Fig there is only one root in the left at low temperatures, which corresponds 
to the sequence "A" in Fig ©-((HI)- When T grows up beyond a specific temperature Tbm, another two 
roots (corresponding to sequences of "B" and "C") will appear. It is noted that the sequence "B" is 
unstable. Tbm increases with Cq and decreases with Ce, as shown in Fig 0. Data here are normalized 
by the Stephen-Boltzmann Limit of ideal QGP, which is defined as 

esB = (16+|n/)^T4; (17) 




Figure 5: Ratio of occupation at different Figure 6: Ratio of occupation at different 
Temperatures, where Cq = 0.8 and Ce = Temperatures, where Ca = 0.8 and Ce = 
by.lQ-^°MeV-^. Data below Tc is just to show 1 x lQ-^°MeV-^. 
the intention of the sequences. 

It is easy to prove that the Free Energy of sequence "C" is smaller than sequence "A". That means 
sequence "A" is a metastable state when T > Tbm like Superheat Liquids. There might be an upper limit 



6 





■ 

A 


L. 


-■-A 
- •- B 
* c 




■ 

: / 

■ 


1 


■ 

\ r 

1 ' 1 ■ 1 ' 1 ' 1 



100 200 300 400 500 600 



T (MeV) 



Figure 7: Energy density of the QGP with 
binary bound states at different Temperatures, 
where C„ = 0.8 and Ce 5 x IQ-^^MeV-^. 
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Figure 8: Energy density of the QGP with 
binary bound states at different Temperatures, 
where = 0.8 and Ce 1 x IQ-^^MeV-^. 



of temperature in the strongly coupled Quark Gluon Plasma to melt out the bound states somewhere at 
T > Tbm. 

According to Eq H13|) and H12|) . one could find (1 — Ca)Ce is major term in Eq H13|l . That makes Tbm 
nearly a function of (1 — Ca)Ce with a small correction less than IMeV. It can be found clearly in Fig 
© and H1U|) . Detailed investigation on the behavior of Tbm implies that linear to [(1 — Ca)Ce]~^ 

approximately in the region of 0.4 x lO^^MeV^ < [(1 - Ca)Ce]"^ < 5 x lO^^MeV^. Calculations with 
larger parameters than 5 x may contain much larger uncertainties and instabilities. 

Details of the transition between sequences depend on the local fluctuations, which are not clear so 
far. Approximately, 



VleIbJ ^ KVlCsbT^ 



^SB X^^^SB 

where Vl is the concerning local volume, and 



The influence of local fluctuations on sequence "A" could be described as. 



°^ TFTJTT-^^ 7^ (18) 



SB - £A £B-Ca£SB T^/'^[£b / ^SB - Co) 

which increases with temperature ( shown in Fig ). It implies that when T — > oo, the system 
should remain at sequence "C" finally. If the local fluctuation is large enough at Tbm, there might be 
a stage of mixed states from Tbm before the whole system transform to sequence "C". In this stage, 
local states could evolve to each other between sequences "A" and "C" by fluctuations. Otherwise, if the 
local fluctuation is too small to drive the sequence "A" to leave the metastable region at Tbm, a sudden 
conversion may occur somewhere at T > Tbm when the local fluctuation grows large enough. 

As shown in Fig Q and ©, the energy density of the whole system may be larger than those of 
heavy quasiparticles or conventional QGP in some regions. It is even larger than the Stephen-Boltzmann 
Limit when T is close to Tc. Since the detailed behavior of color deconfinement is not discussed here, 
the data close to Tc are just qualitative. One of the obvious raisons to enhance the energy density 
is that more degrees of freedom are inserted into the system. Another raison is that more light and 
invariable current quark masses are used in the calculations instead of those ones proportional to T|19|. 
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Ce (10 MeV ) 
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Figure 9: Tbm dependence on parameters of Figure 10: Tbm dependence on parameter of 
and Ce. (1 - C„)Ce. 



Accordingly, data below Tc are not physical. They are only listed in the figures to show the inclination of 
the sequences. When the temperature increases, the ratio of occupations grows up, which could reduce 
the energy density to a factor near Cq relative to Esb, before the state transition takes place at a limit 
temperature T > Tbm- After the conversion, bound states will be melted out and the ratio of occupations 
will decrease to zero. The energy density increases and then the system restores to ideal QGP. Thinking 
that Ca may increase with temperature, (because the reduction of phase space in the " ball" is mainly 
subjected to low momentum regions,) the transition may be more smooth than expected. If Ca grows to 
1 somewhere, bound states will exist for a long period and the transition may vanish at all. 

5. Summary 

Considering the volume of bound states, we provided a phenomenological mechanism to discuss its 
influence on the properties of QGP. There are different solutions for the masses and volumes of bound 
states in different conditions [2^. We constructed an extended Bag Model as a simphfication to describe 
the volume effect qualitatively. Although the details of the model may not be serious, it is still helpful 
for us to understand the statistical behaviors of the bound states with excluded volumes in QGP. 

The contribution of the bound states to the system could be described by the solutions of a self- 
consistent equation (|16|l . From its numerical calculation, we found that the energy density near Tc is 
higher than the one with only quarks and gluons, and the system is nearly occupied by the bound states. 
The ratio of occupations increases with T, and the energy density relative to ess or ideal gas is reduced. 
When the temperature is larger than Tbm, multi-roots of Ea H16|l start to appear. A conversion takes 
place in the system to melt all the bound states at the limit temperature T > Tbm- 
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